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Abstract 



■ Recent developments in studies of multiparticle correlations in high 

■^j- \ energy particle collisions are reviewed. Both experimental data and 

C*") ■ theoretical results in quantum chromodynamics are discussed. Appli- 

cation of the developed methods to the coherent, squeezed and cor- 
related states of photons is considered. Some speculations concerning 
possible coherent and squeezed states of pion fields as well as their 



■ specific features are described. 

p^. 1. INTRODUCTION 

^h. The study of correlations of quantum fields is a common subject in statistical 

r~| ! physics, quantum optics, multiparticle production. Therefore, it is not at all 

surprising to find out common methods of research in those fields. I shall 
describe some newly obtained results in studies of correlations of particles 
produced in high energy collisions as well as their relation to the correlations 
of photons and, in particular, to the squeezing phenomenon. 

The nature of any source of radiation (of photons, gluons or other enti- 
ties) can be studied by analyzing multiplicity distributions, energy spectra, 
various correlation properties, etc. A particular example is provided by the 
coherent states of fields which give rise to Poisson distribution. However, 
in most cases one has to deal with various distributions revealing different 
dynamics. In particle physics dealing with strong interactions, one relies on 
quantum chromodynamics (QCD) as a substitute of quantum electrodynam- 
ics (QED) in case of photons. The fields of quarks and gluons are considered 
instead of electrons and photons. Let us stress that quarks are described by 
spinors similar to electrons, and gluons are massless vector particles like pho- 
tons. The important difference is that gluons carry the color charge while the 
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photons are electrically neutral. It gives rise to the self-interaction of gluons 
as well as to famous properties of the asymptotic freedom and confinement. 

Let us turn now to some experimental facts about high energy particle 
interactions. When two high energy particles collide, the bunch of new parti- 
cles is produced in each event. It is common to plot the distribution of those 
events in the number of particles produced (the multiplicity distribution, for 
short). Usually, they differ drastically from the Poisson distribution, clearly 
indicating strong correlations. Several years ago the experimentalists of UA5 
Collaboration in CERN noticed Q a shoulder in the multiplicity distribu- 
tion of particles produced in pp collisions at energies ranging from 200 to 
900 GeV in the center of mass system. (Let me remind here that the mass of 
the proton is about 1 GeV, i.e. much smaller.) It looked like a small wiggle 
over a smooth curve and was immediately ascribed by theorists to processes 
with larger number of Pomerons exchanged in the traditional schemes. More 
recently, several collaborations studying e + e~ collisions at 91 GeV in CERN 
reported (see, e.g., (2|, H|) that they failed to fit the multiplicity distributions 
of produced particles by smooth curves (the Poisson and Negative Binomial 
distributions were among them). Moreover, subtracting such smooth curves 
from the experimental ones they found steady oscillatory behaviour of the 
difference. It was ascribed to the processes with different number of jets. 

The new sensitive method of theoretical analysis of distributions was 
proposed in [|], [j| (for the review see ||). It appeared as a byproduct of the 
solution of the equations for generating functions of multiplicity distributions 
in quantum chromodynamics (QCD). It appeals to the moments of the mul- 
tiplicity distribution. According to QCD, the so-called cumulant moments 
(or just cumulants), described in more details below, should reveal the os- 
cillations as functions of their ranks, while they are identically equal to zero 
for the Poisson distribution and are steadily decreasing positive functions 
for Negative Binomial distribution so widely used in phenomenological anal- 
ysis. Experimental data show the oscillatory behaviour of cumulants (see 
J?]]) of multiplicity distributions in high energy inelastic processes initiated 
by various particles and nuclei, even though some care should be taken due 
to the high multiplicity cut-off of the data. When applied to the squeezed 
states, the method demonstrates || the oscillations of cumulants in slightly 
squeezed states and, therefore, can be useful for their detection. In such 
circumstances one is tempted to speculate about the alternative explanation 
when considering possible similarity of these findings to typical features of 
squeezed and correlated states. However, first let us describe briefly the 



method of analysis of the distributions which we rely upon. 



2. DISTRIBUTION FUNCTION AND ITS MOMENTS. 

Here I summarize definitions and notations for the values that are used to 
characterize various processes of inelastic scattering according to number of 
particles produced in these processes. The relations connecting these values 
to each other and the examples related to some distributions typical in the 
probability theory are provided. 

Any process of inelastic scattering (that is the scattering with new par- 
ticles produced) can be characterized by the function P n , the multiplicity 
distribution function. The value of P n denotes the probability to observe n 
particles produced in the collision. It is clear that P n must be normalized 
to unity: 

oo 

E P n = l- (2.1) 

n=0 

Sometimes the multiplicity distribution of particles produced can be conve- 
niently described by its moments. It means that the series of numbers P n is 
replaced by another series of numbers according to a certain rule. All these 
moments can be obtained by the differentiation of the so-called generating 
function G(z) defined by the formula: 

oo 

G(z) = p n* n - (2-2) 

ra=0 

Thus, instead of the discrete set of numbers P n we can study the analytical 
function G(z). 

We will use the factorial moments and cumulants defined by the following 
relations: 
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where 
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(n) = J2 Pnn (2.5) 

n=0 



is the average multiplicity. 

Reciprocal formulas expressing the generating function in terms of cu- 
mulants and factorial moments can also be obtained: 
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(F = Fi = 1), (2.6) 
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The probability distribution function itself is related to the generating func- 
tion in the following way: 
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Factorial moments and cumulants are connected to each other by the follow- 
ing recursion relation: 

5-1 

Fg = ^ Cg^Kg^mFm, (2.9) 
m=0 

where 

9-1 m!(g — m — 1)! 

are binomial coefficients. Relation (2.9) gives the opportunity to find the 
factorial moments if cumulants are known, and vice versa. 

It must be pointed out that cumulants are very sensitive to small varia- 
tions of the distribution function and hence can be used to distinguish the 
distributions which otherwise look quite similar. 

Usually, cumulants and factorial moments for the distributions occuring 
in the particle physics are very fast growing with increase of their rank. 
Therefore sometimes it is more convenient to use their ratio: 

H q = ^, (2.10) 



which behaves more quietly with increase of the rank q remaining a sensitive 
measure of the tiny details of the distributions. 

In what follows we imply that the rank of the distribution function mo- 
ment is non-negative integer even though formulas (2.3) and (2.4) can be 
generalized to the non-integer ranks. 

Let us demonstrate two typical examples of the distribution. 

1. Poisson distribution. 

The Poisson distribution has the form: 

( n \n 

P n = { -L-exp(-(n}). (2.11) 

Generating function (2.2) can be easily calculated: 

G(z)=exp({n)z). (2.12) 
According to (2.3) and (2.4) we have for the moments of this distribution: 

F q = 1, K q = H q = 5 ql . (2.13) 

2. Negative binomial distribution. 

This distribution is rather successfully used for fits of main features of 
experimental data in particle physics. It has the form: 

T(n + k) f(n)\ n f (n)\- (n+k) , , 

p " - r ( „ + nrm (x) ( 1 + x) • < 2 - 14 > 

where T is the gamma function, and A; is a fitting parameter. 

At k = 1 we have the usual Bose distribution. The Poisson distribution 
can be obtained from (2.15) in the limit k — > oo. 

Generating function for the negative binomial distribution reads: 

G(,)=fi-^y', 
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and the moments of this distribution are: 

= T(k + q) 
q T{k)ki ' 
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3. QCD EQUATIONS AND MOMENTS OF 
DISTRIBUTIONS. 



The multiparticle production processes are described in quantum chro- 
modynamics as a result of the interaction of quarks and gluons which leads 
to creation of additional quarks and gluons forming the observed hadrons 
at the very last stage. The most typical features of the processes are de- 
termined by the vector nature of gluons and by the dimensionless coupling 
constant. The gluons are colour charged in distinction to photons which 
have no electric charge. Therefore, they can emit gluons in addition to 
quark- antiquark pairs. That is why both quark and gluon jets are consid- 
ered in quantum chromodynamics as main objectives. Their development is 
described by the evolution equations. The main parameter of the evolution 
is the opening angle of the jet or its transverse momentum. The subsequent 
emission of gluons and quarks fills in the internal regions of the previously 
developed angular cones so that they do not overlap (angular ordering) . This 
remarkable property can be exploited to formulate the probabilistic scheme 
for the development of the jet as a whole. Then its evolution equations re- 
mind the well-known classical Markovian equations for the " birth-death" (or 
"mother-daughter") processes. (The detailed discussion of that approach, 
based on the coherence phenomenon, see in ||). 

It is quite natural to start our studies with the simplest case of gluo- 
dynamics. There are no quarks in that case, and interactions of gluons are 
considered only. The system of equations degenerates to the single equation 

G'(y) = f 1 dxK(x) 7 2 [G(y + hxx)G(y + ln(l - x)) - G(y)], (3.1) 
Jo 

where G'(y) = dG/dy, y is the evolution parameter, 

7o = 1 ( 3 - 2 ) 

7T 

as replaces a of electrodynamics and the kernel of the equation is 

K{x) = --(l-x){2-x{l-x)}. (3.3) 

It is the non-linear integro-differential equation with shifted arguments in 
the non-linear term which take into account the energy conservation. It 
can be reduced to the equation for the moments ||] and solved. The most 
prominent feature of the solution are oscillations of cumulants (or ratio H q ) 



; see Q. It has been confirmed by experiment as shown shown in fig. 1. 
It differs from all previously considered phenomenological distributions and, 
from the mathematical point of view, is interesting since it implies that 
QCD deals with non-infinitely divisible distributions in contradistinction to 
conventional ones. In particular, it forbids Poissonian cluster models so 
popular in physics modelling. 

It is interesting to note that it corresponds to rather smooth distribution, 
i.e. there is no obvious one-to-one correspondence between the shapes of 
distributions and behaviour of their cumulants. 

4. APPLICATION TO PHOTONS. 



The above methods can easily be applied to photons ||. Since the gen- 
erating functions of photon distributions for the squeezed and correlated 



states are known [10] one can use the equations of section 2 and get all the 
required moments. Let us consider the most general mixed squeezed state of 
one-mode light. The generating function of the photon number distribution 
was obtained in 11]: 
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The complex parameter (z) is given by relation 



(z) = 2—2 ((g) + i(p)) . (4.4) 



mn = a pp = Tr [gp 2 ^ - (p) 

v qq = Tr (qq 2 ) - (q) 
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mi2 = a pq = ^Tr [g(pq + qp)\ - (p)(q), (4.5) 



T = Tr m = a pp + a qq 

d = det m = a pp a qq - a pq 

, and (Jij are the well known elements of the dispersion matrix. 

The photon distribution function exhibits an oscillatory behaviour if we 
deal with highly squeezed states ( T = a pp + a xx 3> 1 ) for large values 
of the parameter z. A question arises: is it possible to obtain a similar 
"abnormal" behaviour of other characteristics of the photon distribution, 
namely, introduced in section 2 cumulants, factorial moments and their ratio 
H q ? If yes, then in what region of parameters can such anomalies reveal 
themselves? 

The direct differentiation of function In G(u) at it = 1 yields the cumu- 
lants (see section 2) 
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Now let us consider the case of the slightly squeezed state, y= (T-l) <1 
when the photon distribution function does not oscillate. Impose also an 
additional condition 



\z\ 2 



that makes possible to obtain approximate formulas for the functions K q , 
F q , and H q . For K q , we have the following approximate expression: 



K q = g!(_i)9-V-« • (4.7) 

Then recursion relation (2.9) yields: 

F q = q\{-lf 1 - q L-\ 1 ) , (4.8) 

where L~ 1 (x) are generalized Laguerre polynomials. For H q , with g<C7we 
have: 

H q = K q /F q = -— I— » (_i)*+-i g!7 i-9 « 1 . (4.9) 
£<7 (7) 

(If 7 3> <7 the term with highest power of 7 dominates over the rest of sum 
in 1^(7), and F q — > 1 as for Poisson distribution). The exact shape of 
the function iT g is shown in fig.2 a . The distribution function P n does not 
oscillate (fig.2 6 ). 

However, the most abrupt oscillations of the functions K q and H q have 
been obtained when (T — 1) <C 1, but condition 7 3> 1 is not valid. The cor- 
responding curves are shown in figs. 3°, 3 b . Note that the photon distribution 
function is smooth again being approximately equal to zero at q 7^ 1. 

The most regular oscillating patterns of K q and H q are seen at (T—l) ~ 
0.1, \z\ ~ 1 (figs. 4°, 4 b ). Finally we consider the opposite case when the 
photon distribution function P n exhibits strong oscillations while K q and 
H q behave smoothly. Such a behaviour is typical at T ~ 100, \z\ ~ 1 when 
K q exponentially grows while H q monotonically decreases with q (fig. 5). 

Thus we have shown that the cumulants of the photon distribution func- 
tion for one-mode squeezed and correlated light at finite temperature pos- 
sess strongly oscillating behaviour in the region of slight squeezing where the 
photon distribution function itself has no oscillations. And vice versa in the 
region of large squeezing, where the photon distribution function strongly 
oscillates, the cumulants behave smoothly. Hence, the behaviour of cumu- 
lants may provide a very sensitive method of detecting very small squeezing 
and correlation phenomena due to the presence of strong oscillations. 

Let us note also that these methods were successfully applied to predict- 



ing some regularities in behaviour of lasers in the lasing regime [12]. 



5. COHERENT AND SQUEEZED STATES IN PION 

PHYSICS. 



The notion of coherent and squeezed states is not widely applied in 
pion production yet even though several papers deal with the subject [13|, 
[ b§ , |f] , @] , , |H , H , m , H , H . The most striking difference with pho- 
ton case appears due to the requirement of isotopic spin conservation inher- 
ent in pion physics but absent for photons. It gives rise to strong charge 
asymmetry of produced events. If the isospin zero projection of the coherent 
state is considered, it predicts many events without neutral pions, for exam- 
ple. It is demonstrated in fig. 6 borrowed from It used to be related to 
the so-called Centauro-type events observed in cosmic rays. 

Let us consider for definitiveness the process of the collision of two high 
energy nucleons, where in the final state, apart from the two nucleons, a 
certain number of pions has been created. The total isotopic spin of the 
pion system is limited, according to the conservation laws, to the values 
/ = 0,1,2, whereas in the general case it could take values up to n to t, 
where ntot is a total number of pions. This fact significantly affects the pion 
charge distribution. Even if the distribution in the total pion number ntot is 
poissonian (as it happens when pions are produced by classical currents), it 
turns out, that the separate distributions of the charged (n+,ra_) or neutral 
(no) pions are much wider than the Poisson one (see fig. 6). The experimental 
confirmation of this fact could be the "Centauro" events, where the charged 
particles are noticeably dominating, or "anti-Centauro" ones with a large 
number of neutral pions. 

Let us illustrate the idea by considering, following the quasiclassical ap- 
proach of refs. [ 15 1 , [ 18| , the production of many pions in a system with zero 
isotopic spin. The characteristic initial assumption is a possibility of de- 
scribing the pion system that radiates the final state pions as a classical field 
, i.e. that the number of pions per phase space cell is assumed to be big. 
According to the standard reduction formula, the amplitude of generation 
of N pions by the source J equals 

A ai '-' an (k 1 ,...,k n )= lim f Dir a f DJ a W[J] 

N 

exp(iS[7r] + i J d A xn a J a ) n J d A x n e ik " Xn (-dl n - m£)7r°"(x n ), (5.1) 



n=l ' 



where the functional integration over J corresponds to the averaging over 
the characteristics of the pion source. Let us notice, that the radiation of 



a classical current exactly reproduces the language of coherent states. The 
specific feature of pion fields is that one projects the coherent states onto 
the states with a definite isotopic spin, and it gives rise to drastical change 
of final charge distributions. The quasiclassical estimate of the amplitude, 
performed in the assumption on the axial symmetry of the initial interaction 
and on the isotopic symmetry of the pion system (i.e. of the zero total 
isospin) fl4T| , fll5| , fl8|l , leads to the two characteristic conclusions that also 
appear in other publications on this topic. 

Firstly, only the distribution over the total number of pions is poissonian. 
The distributions over the number of neutral and charged pions are much 
wider (see fig. 6). For the state with zero isospin a probability of finding 2n 
neutral pions in the system of 2N pions has a form [13], [19]: 



At large n, N we have a characteristic distribution |23| : 

P(n,N) ~ (n/N)- 1/2 . (5.3) 

Secondly, the conservation of the total isospin leads, for example, to the 
specific angular correlations between the particles with different charges. 
Let us give the characteristic formula for the correlation over the azimuthal 
angle <p, obtained in [24] for the pions having zero rapidity: 



/ do 7 ' ; 9 da 7 * da 71 \ „ 
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— COS (<fi - <fi 2 ). (5.4) 



dki dk 2 



The experimental verification of such predictions is to our opinion very in- 
teresting. 

The similar procedure of projecting squeezed states onto the states with 
definite isospin has been attempted in p5] ]. It produced some interesting 
modifications of fig. 6 which we have no space to discuss here. 

In the recent literature a "disoriented chiral condensate" was widely dis- 
cussed as a possible asymmetry source in the production of charged and 
neutral pions in some fraction of the events ( in particular, in the above- 
mentioned " Centauros" ) . Let us remind, that the transformation properties 
of mesons with respect to the chiral transformations are determined accord- 
ing to the corresponding properties of the order parameter, which character- 
izes the spontaneous breaking of chiral symmetry and is given, for example, 



by the average from the bilinear combination of quark fields 



$ ~ {QLqR), (5.5) 

where qRUS are the right (left) states of the massless quarks. For investi- 
gating the character of the singularity of thermodynamical functions in the 
vicinity of the phase transitions, it is desirable to find a solvable model hav- 
ing the same symmetry. Then, according to the universality principle based 
on the scale invariance near the critical point, the solutions of this model will 
have the same set of singularities. In such an approach the order parameter 
( |5.5[ ) can be rewritten in terms of a set of hadron fields having the same 
symmetry. Therefore, the multipion states become related to the massless 
quark fields and quark-gluon plasma, i.e. there appears a hope to describe a 
phase transition between these so differing phases within such a model. In 
the realistic case of two massless quark flavors the chiral field can be written 
in the form 

$~cr«l-Hf«7T, (5.6) 

where a, tt are the real fields, f are the standard Pauli matrices, the ir - 
meson fields tt form an isotriplet and a is an isosinglet. 

In the case of SU(3) - algebra the number of such fields is already 18. 
They form the scalar and pseudoscalar nonets. 

The dynamics of these degrees of freedom can be described by the la- 
grangian of the linear a model 

L= l -{daf+ l -(d^f-V{a^l (5.7) 

where V is a potential depending on the combination cr 2 +7f 2 . In the standard 
version [^] the spontaneous symmetry breaking occurs via the formation of 
the nonzero vacuum average of the field a. The isotriplet fields remain 
massless, i.e. the pions are the goldstones of the chiral group. 



Let us now assume [17], that in some region of space the vacuum orien- 



tation is different from the standard one, and, for example, 

W) = ZttCosA, (tt) = f n nsm6, (5.8) 

where f n = 93 MeV, and n is a unit orientation vector of tt. Such an 
assumption presupposes a specific scenario of the process, which is still not 
studied in details. 



If the field <I> is isotropic with respect to a direction on the 3-dimensional 
sphere in the 4-dimensional space with the angles defined as 



(<t, 7T3 , 7T2 , 7i"i ) = (cos 6, sin 6 cos (f>, sin 9 sin <f> sin 77, sin 6 sin cos 7?) , (5.9) 
then for the probability distribution of a given state r = cos 2 </> we have: 

o 1/2 

/ drP(r) = — dr] d6 sin 2 6» / #sin0, (5.10) 

Jri JO JO JaxccosrJ 



as it was obtained for the first time in refs.|15|, 23 1: 



P(r) = ^, (5.11) 



i.e. we have again returned to the formula ( |5.2| ). Thus the probability of 
finding an event with a fraction of the neutral pions being smaller than a 
certain given tq is 

W(r < r ) = y/fZ, (5.12) 

which constitutes 10 % even at ro = 0.01. The charge fluctuations in such 
a system are much larger than those that would follow from the Poisson 
distribution, where the distributions are concentrated around r=l/3. 

Let us note that the model of creation of N isoscalar pion pairs is de- 
scribed as 

|*) oc (2a+a± - (4f) N \0), (5.13) 

which is reminiscent of squeezed states projected onto zero isospin but is 
more restrictive than that projection. (Here a+ _ are the creation operators 
of positive, negative, neutral pions, correspondingly). 

6. CONCLUSIONS 



We have shown that the correlation analysis can be done by similar 
methods for photons and pions. Some further understanding of the relation 
between distribution behaviour and features of its moments is awaited for. 
More clear insight into usefulness of the idea of coherent and squeezed states 
as applied to particle physics is required. However, the very first attempts 
in these directions described above are very encouraging. 



References 

[1] UA5 coll., Alner G.J. et al., Phys. Rep. 154 (1987) 247. 

[2] DELPHI coll., Abreu P. et al., Z. Phys. C50 (1991) 185. 

[3] OPAL coll., Acton PD. et al., Z. Phys. C53 (1992) 539. 

[4] Dremin I.M., Mod. Phys. Lett. A 8 (1993) 2747. 

[5] Dremin I.M., Phys. Lett. B 313 (1993) 209; 

Dremin I.M. and Hwa R.C., Phys. Rev. D 49 (1994) 5805. 

[6] Dremin I.M., Usp. Fiz. Nauk 164 (1994) 785; Physics -Uspekhi, 37 
(1994) 715. 

[7] Dremin I.M., Arena V., Boca G., et al., Phys. Lett. B 336 (1994) 119. 

[8] Dodonov V.V., Dremin I.M., Polynkin P.G., and Man'ko V.I., Phys. 
Lett. A 193 (1994) 209. 

[9] Dokshitzer Yu.L., Khoze V.A., Mueller A.H. and Troyan S.I., Basics of 
perturbative QCD, Gif-sur-Yvette, Editions Frontieres, 1991. 

[10] Dod onov V.V., Drem in I.M., Polynkin P.G., Man'ko O.V., and Man'ko 
V.I., [hep-ph 9502394 . 

[11] Dodonov V.V., Man'ko O.V., and Man'ko V.I., Phys. Rev. A 49 (1994) 
2993. 

[12] Hwa R.C., Phys. Rev. C 50 (1994) 383. 

[13] Horn D., Silver R., Ann. Phys. (N.Y.) 66 (1971) 509. 

[14] Botke J.C., Scalapino D.J., Sugar R.L., Phys. Rev. D9 (1974) 813; D10 
(1974) 1604. 

[15] Andreev I.V., Pisma v ZhETF 33 (1981) 384 (JETP Lett. 33 (1981) 
367). 

[16] Vourdas A., and Weiner R.M., Phys. Rev. A 36 (1987) 5866. 

[17] Bjorken J.D., Acta Physica Polonica B23 (1992) 561. 

Bjorkcn J.D., Kowalski K.L., Taylor C.C., SLAC-PUB-6109 (1993). 



[18] Anselm A.A., Phys. Lett. B217 (1989) 169. 

[19] Kowalski K.L., Taylor C.C., CWRUTH-92-6, [hep-ph/921 12821 (1992). 

[20] Rajagopal K., Wilczek F., Nucl. Phys. B379 (1993) 395; ibid B404 
(1993) 577. 

[21] Kogan LI., Phys. Rev. D48 (1993) 3971. 

[22] Gavin S., Mtiller B., Phys. Lett. B329 (1994) 486. 

[23] Anselm A.A., Ryskin M.G., Phys. Lett. B266 (1991) 482. 

[24] Blaizot J.P, Krzywicki A., Phys. Rev. D46 (1992) 246. 

[25] Dremin I.M., and Hwa R.C. (in preparation). 

[26] Gell-Mann M., Levy M., Nuovo Cimento 16 (1960) 705. 



